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Abstract 



By examining two counterexamples to the existing theory, it is 
shown, with mathematical rigor, that as far as scattered particles 
are concerned the true distribution function is in principle not deter- 
minable (indeterminacy principle or uncertainty principle) while the 
average distribution function over each predetermined finite velocity 
solid-angle element can be calculated. 

1 Introduction 

The distribution function, playing the most fundamental and most central 
role in classical statistical mechanics, is defined by the limiting process: 



where Ar = AxAyAz is a position volume element, Av = Av x Av y Av z a 
velocity volume element and AN the number of the particles found in the 
phase- volume element ArAv. This concept has served us since the start of 
statistical mechanics and no serious challenge was put forward against it. 

However, examining the definition ([T]) with mathematic rigor leads us 
to stimulating ideas. For instance, if, for whatever reason, some dimensions 
of ArAv have to be kept finite (not infinitesimal), we have no choice but 
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to regard them as the uncertainty ranges over which f(r,v,t) expressed 
by (H|) is not truly determinable. To exclude such uncertainty, we need to 
show that the ratio in (pQ) will tend to a definite value no matter in what 
way Ar and Av approach zero. This is by no means an insignificant or 
dispensable task in view of the fact that many intriguing examples in math 
have manifested that a multi-variable function should be considered to have 
no limit wherever it has path-dependent limits [lj. (Also note that taking 
partial derivatives at such places becomes unjustified and improper.) 

Being exposed to a variety of thought experiments and doing thorough 
examinations on them [21 [3], we are now convinced that there indeed exist 
inherent constraints to limit how Ar and Av tend to zero; or, in other words, 
there indeed exists a certain type of uncertainty principle. The situation, to 
some extent, resembles what happens to the uncertainty principle of quan- 
tum mechanics, which was in history, still is in usual textbooks, obtained 
from analyzing experimental facts. 

Since the studies of distribution function were strongly influenced by 
Boltzmann's initial approach on the Boltzmann gas (by the term it is meant 
that particles there interact with each other via binary short-range forces), 
we shall here concern ourselves with the Boltzmann gas only. 

In section 2, two plain counterexamples to the existing theory are pre- 
sented, which hints that the commonly held notion about distribution func- 
tions is actually unsound. In section 3, the two counterexamples are formu- 
lated with help of the textbook methodology of collisions; it is found that 
only the average distribution function over each finite velocity solid-angle 
element can be calculated, while the true distribution function is in principle 
not determinable. In section 4, a special type of indeterminacy principle, or 
uncertainty principle, is addressed and discussed. Section 5 concludes this 
paper. 

2 Counterexamples to the existing theory 

To hint that the conventional notion about distribution functions indeed 
involves difficulty, two plain counterexamples against the Boltzmann theory 
are addressed in this section. 

The first counterexample is related to the setup shown in Fig. 1 , in which 
there are two dilute particle beams, produced by two sources, colliding with 
each other. On the premise that the initial distribution functions of beam 
1 and beam 2 are respectively given as 

/((r,vi,t) and ^(r,v^,t), (2) 
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we try to determine the distribution function of the beam 1 particles scat- 
tered by beam 2 particles, denoted as f(r, v\,t). In the consideration herein, 
/{ and j'i are assumed to be quite general, containing no <5-functions to 
avoid controversy [5], which can be accomplished by letting the two particle 
sources have finite temperatures and finite outlets. For / = f(r, V\,t), we 
write Boltzmann's equation in the form[5]: 
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where m is the mass of a beam 1 particle, F the external force, and (Sf/St) gSu i n 
and {Sf /St)\ oss represent the collision terms making / increase and decrease 
respectively. 



Collision region 



fl 



r t i 
• r 



Fig. 1: Two particle beams colliding with each other. 

Inside the collision region, due to the aforementioned assumption that 
the two beams are dilute, / must be much smaller than /{ or f' 2 and thereby 
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which means that secondary collisions can be neglected. By adopting that 
the situation is stationary (the sources produce particles in a constant man- 
ner during the time of interest) and there is no external force acting upon 
the particles, we arrive at 
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Outside the collision region, we get 
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In Fig. 2, what equations © and © say about the distribution function 
versus the distance to the system's center is schematically depicted by the 
dotted curve. 
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Fig. 2: Schematic of the distribution function versus the 
distance to the system's center. The dotted curve is related 
to Boltzmann's equation while the solid one to the reality. 

However, our physical intuition, as well as any well-conducted exper- 
iment, tells us something remarkably different: the distribution function 
will, beyond a certain distance to the system's center, diminish gradually to 
zero, as the solid curve in Fig. 2 illustrates. 

The second counterexample is shown in Fig. 3, in which the particles of 
a dilute beam (of (5-function type or not) collide with a solid boundary. 



• r 

Fig. 3: A particle beam colliding with a solid-surface bound- 
ary 

Again, Boltzmann's equation is employed to study the distribution func- 
tion of the particles scattered by the boundary, denoted as /(r, v±,t) = f. 
If the situation is stationary and no external force applies, the distribution 
function obeys 




T 
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In Fig. 4, what equation ([7]) says is given by the dotted curve while what 
our physical intuition, as well as real experiments, says by the solid curve. 



»- r 

Fig. 4: Schematic of the distribution function versus the 
distance to the boundary. The dotted curve is related to 
Boltzmann's equation while the solid one to the reality. 

In view of the fact that all moving particles can be deemed as ones 
scattered either by other particles or by boundaries, the foregoing two coun- 
terexamples, though somewhat heretical, possess general significance. 

3 Formulations of the counterexamples 

In order to understand the counterexamples given in the last section, we 
shall try to formulate the distribution functions by all means available to 
us. 

r' 

Detector I 1 /U^A^! 
(a) 

Fig. 5: (a) The detector's opening is placed toward the colli- 
sion region; (b) when the opening shrinks to zero, a clear-cut 
effective region is formed. 

For the first counterexample, referring to Fig. 5, a virtual detector is 
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placed at the point r where the velocity distribution is of interest and the 
opening of the detector is toward the collision region. We then assume 
that every beam 1 particle entering the detector within the velocity element 
v^AwiAili (y\ = |vi| and Q\ is the solid-angle of vi in the velocity space) 
will be registered as an "effective particle" . According to the standard the- 
ory, the distribution function is the limit of the following ratio: 



in which Ar is a spatial volume element just inside the detector opening and 
AiV"i is the number of all the effective particles in Ar. It turns out that ([8]) 
can be directly and readily calculated except that there is a compromise. 

The compromise is related to the fact that © can be evaluated only 
under the assumptions: 



In what follows, we shall do the calculation with Q adopted, and later on 
investigate what happens if Af2i also shrinks to zero. 

Let ASo be the area of the detector opening and AUq be the solid angle 
domain formed by A So and a point r' in the collision region, as shown in 
Fig. 5(a). In order for the detector to detect the ratio (jHJ) under ([9]), ASo, 
and thus Af2o> needs to be infinitely small. With AOi being finite and 
AOo being infinitesimal, it is found that the collisions occurring inside the 
region enclosed by the spatial cone — A^i will possibly produce effective 
particles while the collisions taking place elsewhere will not, as shown in 
Fig. 5(b). For convenience, the region, defined by r and — AOi, is named as 
the effective region and the path linking from a point in it to the detector 
opening as an effective path. 

Following the textbook methodology of treating collisions almost exactly [5j, 
the number of all relevant effective particles can be represented by 



At' / dr> / dv[ / dv' 2 / dn ■ 2ua(n)f[(r'y i ,t')ft ! (r'y 2 ,t'), (10) 



where / dr' ■ ■ ■ runs over the entire effective region (over all effective paths 
in the sense), 2u = \v[ — v' 2 \, cr is the cross section, J d£l ■ • • runs over all the 
passages through which effective particles move and enter the detector, and 
t' = t — |r — r'\/v\ reflects the delay between the collision time at r' and the 
arrival time at r. 

Since cx(Q) and dQ, are conveniently defined in the center-of-mass frame 
while the passages from the collision location to the detector opening are 



AN, 



(8) 
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Ar -> 0, Avi -> and AOx 0. 
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conveniently denned in the laboratory frame, we must do something other 
than what textbooks have elaborated. With the aid of the notations c' = 
(vi + v' 2 )/2, c = (vi + v 2 )/2, u' = (vi - v' 2 )/2, u = (vi - v 2 )/2 and 
u = |u'|, the conservation laws of energy and momentum can be written as 
(every particle of beam 1 and beam 2 is assumed to have the same mass) 

c' = c and u = |u | = |u|. (11) 

By virtue of u = |u'|, expression (|lUp becomes 

At' J dr' J dc' j dn' J u 2 du J dfl- \\J\\2ua(n)f[f^ (12) 

where ||J|| = ||d(vi, v' 2 )/d(c' , u')|| = 8, / d£l' ■ ■ ■ runs over all possible direc- 
tions of v[ — v 2 , /{ and / 2 stand for, respectively, 

/{(r', c' + u',0 and ti(r',c'-u',t'). (13) 



c 



Vl 




Fig. 6: Relation between c, u and vi in the velocity space. 

The geometry in the velocity space, shown in Fig. 5(a) and Fig. 6, informs 
us that 

J u 2 du J dQ--- = vfAv^fio (14) 

Fig. 5(a) also shows that the particles will, at the detector opening, occupy 
the spatial volume: 

Ar = r - r'\ 2 An v 1 At'. (15) 
Taking (fT2|) . (fTl"]) and (fT5|) into account, we can express ©, under ([9]), as 

f(r,v u AIM) = dr'J dc'J d^-^§^. (16) 

All quantities in (|16p . explicit or implicit, are well defined. For instance, 
u = (fino — c') with no = (r — r')/\r — r'|, |u'| = [u| = u while u' points in 
the direction of 0', and f2 is defined by u and u'. 
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It is obvious that expression (|16|) represents the distribution function 
averaged over the given, or chosen, finite velocity solid-angle element A$7i. 
At this stage, several points should be made, (i) No approximation has been 
introduced in the formulation, (ii) This formula is consistent with the solid 
curve in Fig. 2; interested readers may confirm it in real or computational 
experiments, (iii) The formalism lives with discontinuity comfortably, (iv) 
If necessary, the secondary collisions along each effective path can be taken 
into account in a relatively easy manner [3]. 

We now look at whether or not 'the true distribution function' can be 
determined. As mentioned in the introduction, if the limit of ([8]) exists, we 
shall meet with it as Ar, Av± and Af^i tend to zero in whatever way. That 
is to say, the limit of (|16p . namely 



expresses nothing but the true distribution function. To see whether (|17p is 
mathematically meaningful, we inspect our formulation from ([8]) to (|17p and 
examine each of the involved limiting processes. The following three groups 
of limiting processes appear very essential, (i) ASq — > 0, thus Af^o — > 0, 
that arises from Ar — > 0. (ii) AOi — > that arises from Avi — > 0. (iii) The 
integration itself can be deemed as a limit: the effective region is divided into 
many small elements in the analysis and the final integration is performed 
continuously over the region. It turns out that these limiting processes do 
not get along very well. As one example, if AVti becomes, when shrinking 
to zero, comparable to AQ.Q, the very definitions of the effective region and 
effective path will break down and the whole formalism will collapse. As 
another example, if we let Af2i approach zero much slower than Af^o does, 
the definition of the effective region manages to hold and expression fjlTf) will, 
with help of dr' ~ |r — r'\ 2 AQidl, become a line integral which corresponds 
to Boltzmann's equation quite well. This outcome, seemingly desirable and 
desired, is actually misleading. Not only that the very spirit of the true 
distribution function, expressed by ([I]), prohibits any ordering between Ar 
and Avi (one type of path-dependence), but that the line integral yields, in 
terms of dealing with this example, a curve that is almost the same as the 
dotted one in Fig. 2. 

We now turn our attention to the second counterexample in the last 
section. Interestingly, we shall encounter roughly the same difficulty. 

Based on the notion that particles scattered by boundaries behave very 
much like ones emitted by boundaries [6], we describe the scattered particles 




(17) 
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of this setup in terms of the emission rate p, which is defined as such that 



dN = p(t', r', v', n')dt'dS'v' 2 dv'dn' (18) 

represents the number of the particles scattered in dt' from a boundary 
patch dS' (located at r') and emerging within the velocity range v l2 dv'dVL' . 
Generally speaking, p depends on many factors, such as the distribution 
function of the incident particles, surface's geometry and surface's property, 
and p should be determined with help of experiments. For the purpose of 
this paper, we shall do our formulation on the premise that p is completely 
known. 

As shown in Fig. 7, a virtual detector is placed at r; ASo and Af^o are 
defined in the same way as in the first counterexample. Furthermore, the 
particles entering the detector within the velocity element v\ Av\ASl\ are 
also registered as effective particles and the assumptions expressed by Q 
are also adopted at the beginning. 



AS" 




r -JVU-ASb 



Detector 



Afii 



Fig. 7: A virtual detector is placed toward the boundary. 

The area on the boundary surface enclosed by — Aili (starting from r) 
will be called the effective area, denoted as AS'. The particles emitted from 
a point r' on AS' will, at the detector opening, occupy the spatial volume 
|r — y'\ 2 AVL^v' dt 1 . Integrating all the contributions from the effective area, 
we obtain "the distribution function" as 



p(t',r',v',n')dt'dS'v' 2 dv'dn' 
AS',Aui,An |r — r'\ 2 Aflov'dt' ■ vjAviAfli 



With help of AS -» (thus A°, -> 0) and Avi -> 0, expression (fTU|) 
becomes 

fU AO ^ 1 f p(t',r',v l7 n')dS' 
v 1 AUi J as' |r - r'\ z 

where t' = t — \r — r'\/vi. 
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Expression (|20p is nothing but the distribution function averaged over 
the finite velocity solid-angle domain. Can we arrive at the exact distribution 
function along the line? The answer to it is "no" again. If Af2i goes smaller 
and becomes comparable to AQq, the effective area will no longer be effective 
and the whole formalism will no longer hold. 

Before leaving this section, it is relevant and interesting to examine 
whether or not we can, for our two foregoing cases, determine the limit: 

ANi 

A/-»o,Au7 : ^o,Ani-*o ASqAI ■ i^AijiAOi ' 

in which A So, representing the opening area of the virtual detector, is kept 
to be finite (while AJ7i — > 0). Expression (f2T|) . though conjugate to (fT6|) . is 
not calculable in light of the fact that, to do the evaluation, we have to divide 
ASo into many dSo and then determine the true distribution function on 
each infinitesimal dSo, which cannot be done with Afii — > 0, as just revealed. 



n AaAt ' An , (21) 



4 A special type of uncertainty principle 

The last section has shown that there are intrinsic and unsurmountable dif- 
ficulties to prevent us from formulating the genuine distribution function 
of scattered particles (unless there exists sure uniformness in space and in 
time); however, by partitioning the velocity solid-angle space into many fi- 
nite elements, the average distribution function over each of the elements can 
be calculated. For a reasonably complex system, the taking-average strategy 
may be thought of as a 'better-than-none' one: once the way of partition is 
predetermined, how the genuine distribution function varies within each of 
the finite velocity solid-angle elements will be indeterminable and all possi- 
ble effects caused by the variation within each of them will be beyond the 
investigation, which literally means that the link from any limited number 
of investigations (one way of partition defines one investigation) to the real 
behavior of the system is not truly certain. The aforementioned concep- 
tion, concerning in what sense we can and cannot determine distribution 
functions, might be regarded as a special type of indeterminacy principle or 
uncertainty principle. 

It is obvious that the integral formalism given in the last section plays 
a vital role in this issue. Not only that it serves as a mathematical demon- 
stration on why the uncertainty principle has to be introduced, but that it 
offers an effective methodology to formulate the Boltzmann gas under the 
uncertainty principle. The following perspectives may help us to compare 
between the standard approach and this approach. 



11 



1. The standard theory is based on the tacit assumption that all distri- 
bution functions are either differentiable or can be approximated by 
differentiable functions, which is, however, in conflict with the fact 
that a variety of real distribution functions are discontinuous or quasi- 
discontinuous (of nonuniform continuity). This approach, while deny- 
ing the legitimacy of differentiating distribution functions, lives with 
discontinuity comfortably: the final formula is in an integral form 
capable of handling discontinuity; and it yields a stepwise varying dis- 
continuous distribution function everywhere. 

2. As well known, one of the assumptions of Boltzmann's equation is the 
time-reversibility of particle collision, while one of the final conclusions 
of Botlzmann's equation is the time-irreversibility of particle collisions; 
the logic paradox there has inspired and continues to inspire a variety 
of discussions in the literature. Whereas, this approach is inherently 
of time-irreversibility in view of the fact that the resultant distribu- 
tion function is obtained by integrating the distribution functions at 
previous times over a region (or area) , via which a lot of historical and 
detailed information is erased explicitly. 

3. In the standard approach, after a distribution function is obtained by 
solving Boltzmann's equation, presumably with help of specific initial 
and boundary conditions, the task of the kinetic theory is finished in 
a once-for-all manner. In this approach, we may alter the partition 
of finite velocity solid-angle elements and restart the calculation. By 
repeatedly doing so, a number of resultant distribution functions for 
the system can be obtained. If the system is sufficiently complex and 
the evolution time is sufficiently long, it is quite possible to find that 
same initial and boundary conditions result in dramatically different 
outcomes. This is partially related to the uncertainty nature of this 
approach and partially related to the non-uniqueness nature of our 
world. 

4. As one of the most important features, this approach produces, in 
many simple cases, definite results that are directly verifiable with 
help of experimental or computational means. In particular, the two 
counterexamples against the standard theory presented in Sect. 2 have 
been nicely treated in Sect. 3. 

Many issues, closely related and not so closely related, have yet to be 
explored; some of them are discussed by the author elsewhere [7] . 
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5 Conclusion 



It is often said that to decide whether a theory bears physical truth, one 
must be able to falsify it. Unfortunately, Boltzmann's theory has been 
established for so long that people in this community lose interest in doing 
so. This paper should, at least, serve as a clue that a considerable number 
of concepts and methodologies of the existing statistical theory need serious 
reconsideration. 
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